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The pseudospin symmetry is a relativistic dynamical symmetry connected with the small com- 
ponent of the Dirac spinor. The origin of pseudospin symmetry in single particle bound states in 
atomic nuclei has been revealed and studied extensively. By examining the zeros of Jost functions 
corresponding to the small components of Dirac wave functions and phase shifts of continuum states, 
we show that the pseudospin symmetry in single particle resonant states in nuclei is conserved when 
the attractive scalar and repulsive vector potentials have the same magnitude but opposite sign. 
The exact conservation and the breaking of pseudospin symmetry are illustrated for single particle 
resonances in spherical square-well and Woods-Saxon potentials. 
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The concept of pseudospin (PS) is often introduced to 
reveal the dynamical nature of quantum systems. More 
than 40 years ago the pseudospin symmetry (PSS) in 
nuclear single particle states was observed: PS dou- 
blets with quantum numbers (n r ,l,j = 1 + 1/2) and 
(n r — 1,1 + 2,j = I + 3/2) are nearly degenerate [1, 2]. 
Since then much efforts had been devoted to explore 
the origin of the PSS (see, e.g., [3-5]) until it was 
shown that the PSS in nuclei is a relativistic symme- 
try which is exactly conserved when the scalar and vec- 
tor potentials have the same size but opposite sign, i.e., 
£(r) = S(r) + V(r) =0 [6]. However, this condition is 
never met in finite nuclei, because in this limit there are 
no bound nuclei any more. Later it was found that the 
PSS is exact under a less strict condition, dS(r)/dr = 0, 
and to what extent the PSS is conserved is related to 
the competition between the centrifugal barrier and the 
PS orbital potential [7, 8]. This condition can also not 
be met in realistic nuclei, therefore experimentally one 
always finds that the PSS is broken. The above men- 
tioned two conditions also result in the spin symmetry 
(SS) in anti-nucleon spectra [9, 10] which is much better 
developed than the PSS in nuclear single particle spec- 
tra [10, 11]. The SS and PSS have been studied exten- 
sively within the relativistic framework in which static 
mean fields dominate, including the PSS in deformed nu- 
clei [12, 13] and the SS for A spectra in hyper-nuclei [14]. 
The relevance of the PSS in nucleon-nucleus and nucleon- 
nucleon scatterings has also been discussed [15-18]. The 
readers are referred to Ref. [19] for a review and Refs. [20- 
25] for some recent progresses. 

In recent years, there has been an increasing interest in 
the exploration of nuclear single-particle resonant states 
especially in the study of exotic nuclei with unusual N/Z 
ratios [26-32]. In these nuclei, the neutron (or proton) 
Fermi surface is close to the particle continuum, thus 
the contribution of the continuum and/or resonances is 
important [33, 34]. The study of symmetries in reso- 



nant states is certainly an interesting topic. There have 
been some investigations of the PSS in single particle 
resonances [35-38]. For example, the PSS for the reso- 
nant states in 208 Pb is investigated by solving the Dirac 
equation with Woods-Saxon-likc vector and scalar poten- 
tials in combination with an analytic continuation in the 
coupling-constant method [26, 27, 39] and it was found 
that the difiuseness of the potentials plays an important 
role in the splitting of energy and width of resonant PS 
partners [35]. However, in all these studies, the PSS in 
resonances was investigated numerically and a rigorous 
justification of the PSS in single particle resonant states, 
like that for bound states given in Ref. [6], is still ab- 
sent. In this Letter, we show that the PSS in single par- 
ticle resonant states in nuclei is also exactly conserved 
under the same condition for the PSS in bound states, 
i.e., E(r) = or dT,(r)/dr = 0. We will also illustrate 
the exact conservation and the breaking of PSS in single 
particle resonances in spherical square-well and Woods- 
Saxon potentials. 

In a relativistic description, nuclei are characterized 
by a strong attractive scalar potential S(r) and a strong 
repulsive vector potential V(r) [33, 34]. The Dirac equa- 
tion for a nuclcon reads 

[a-p + P (Af + S(r)) + V(r)] i/)(r) = e^(r), (1) 

where a and /3 are the Dirac matrices and M is the nu- 
cleon mass. For a spherical nucleus, the Dirac spinor 
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where Yj m 



, y)) is the spin spherical harmonic. F nK (r)/r 
and Gn K (r)/r are the radial wave functions for the upper 
and lower components with n and n radial nodes, k = 
(-l)i+*+V2(j + i/2) andf= Z-sign(re). The radial Dirac 



2 



equation is then derived as 




(F{r) 
\G(r) 



(3) 

where £(r) is defined earlier, A(r) = V(r) — S(r), and e is 
the eigenenergy. For brevity we omit the subscripts from 
F(r) and G(r) whenever no confusion arises. This first 
order coupled equation can be rewritten as two decoupled 
second order differential ones. Here we only write down 
the one for the small component to which the PSS is 
directly connected, 



dr 2 
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M- (r) dr 
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dr 



dr 



M + (r)M-(r) 



G(r)=0, (4) 



where M + (r) = M+e-A(r) and M_(r) = M-e + E(r). 
Equation (4) is fully equivalent to Eq. (3). 

For the continuum in the Fermi sea, i.e., e > M, there 
exist two independent solutions for Eq. (4). (Note that 
the following discussions are also valid for the continuum 
in the Dirac sea.) The physically acceptable solution is 
the one that vanishes at the origin. As usual we define 
the regular solution G(r) as the one that behaves like 
jj(pr) as r -» [40] , 



lim G{r)/jr(jpr) = 1, p = \J t 2 - M 2 . 



(5) 



We now turn to the asymptotic behavior of the reg- 
ular solution as r —> oo. At large r the potentials for 
neutrons vanish and the wave functions oscillate, Eq. (4) 
becomes a Ricatti-Bessel equation with angular momen- 
tum / and the solution can be written as a combination 
of the Ricatti-Hankcl functions, 
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where (p) is the Jost function for the small component 
and h^(pr) the Ricatti-Hankel functions. 

In nuclei the potentials V(r) and S(r) share some gen- 
eral properties, e.g., they are analytic functions of r, van- 
ish when v — y oo . and have no singularities. Under such 
conditions, the Jost function is an analytic function of p 
and can be analytically continued to a large area in the 
complex p plane. Here the structure of the p Riemann 
surface on which the Jost functions are defined is more 
complex than the non-relativistic case. For example, the 
square root in the relativistic energy-momentum relation 
e 2 = p 2 + M 2 creates branching points at p = ±iM, thus 
the corresponding Riemann surface is at least two folds. 

In Fig. I the zeros of the Jost function J^?(p) on the 
complex momentum plane are schematically shown. For 
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FIG. 1. (Color online) Schematic picture of the zeros of the 
Jost function on the complex momentum plane. A cut 
is made on the imaginary axis, from p — iM to infinity and 
back to p = —iM. 



simplicity we show only the first sheet with Re(e)>0 
which contains positive energy bound states and resonant 
states, while the other sheet with Re(e)<0 can be used 
to investigate negative energy ones. These two sheets are 
connected by a cut on the imaginary axis, from p = iM 
to infinity and back to p — —iM. Restricted to the 
first sheet and not too large \p\, J^{p) is a single val- 
ued analytic function of p. The zeros of J^{p) are de- 
noted by full circles (bound states), open circles (reso- 
nant states), and crosses (other zeros), respectively. The 
zeros on the positive imaginary axis of the p plane rep- 
resent bound states of the original eigenvalue problem, 
while the zeros on the lower p plane and near the real 
axis correspond to resonant states. The resonance en- 
ergy E rca and width T I CS are de termined by the relation 
E = E Tes — iT Tes /2 = \J p 2 + M 2 . By examining the zeros 
of the Jost function we can study the bound and resonant 
states on the same footing and many known properties of 
bound states can be generalized to resonances straight- 
forwardly. 

In the PSS limit, Eq. (4) is reduced as 



d2 1(1 + 1) , „ r x „ , s 
—5 - + e - M M+(r) 

dr r z 



G(r) = 0. (7) 



For bound states it is an eigenequation that determines 
the eigenenergy e. While for the continuum e can be any 
value > M and we mainly focus on wave functions and 
their asymptotic behavior. For PS doublets with differ- 
ent quantum numbers k and k' with n' = — k + 1, the 
small components satisfy the same equation because they 
have the same pseudo-orbital angular momentum I [6]. 
In particular, for continuum states, we have G K (e,r) = 
G K r(e,r) for any energy e. Because the definition of the 
Jost function J® (p) only depends on the asymptotic be- 
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havior of the small component, we have J3(p) = J^(p) 
on the positive real axis. This equivalence can be gener- 
alized into the complex p plane due to the uniqueness of 
the analytic continuation. Thus the zeros are the same 
for J^{p) and J^(p) : If there exists a resonant state 
with energy E res and width T ICS and the quantum num- 
ber k, there must be another one with the same energy 
and width and quantum number That is to say, the 
PSS in single particle resonant states in nuclei is exactly 
conserved when the attractive scalar and repulsive vector 
potentials have the same magnitude but opposite sign. 
Certainly if we focus on the zeros of the Jost functions 
of PS doublets on the positive imaginary axis of the p 
plane, we come to the well-known PSS for bound states. 

In scattering theories, one can also determine reso- 
nance parameters from the change of cross section or 
phase shift which give us more insights into the resonant 
phenomena. Next we discuss the PSS in resonant states 
by examining the phase shift. Using the asymptotic be- 
havior of the Ricatti-Bessel functions, one obtains from 
Eq. (7), 



G K (r) cx sin ( pr - y + <5« (jp) 



(8) 



where the phase shift <5J? (p) is related to the Jost function 
through jG(p) = \J^(p)\e~ lS - {p) . Whenever 6%(p) = 
mr + 77 , there is a resonant state and its width is deter- 
mined by the tangent of the phase shift function 8^{p). 
In the PSS limit, the coincidence between G K (r) and 
G K >(r) means that S^(p) = S^,(p) for any value of p. 
Therefore resonance parameters of PS doublets are the 
same. 

It has been postulated in previous numerical studies 
that widths of PS doublets should be different even in 
the PSS limit because centrifugal barriers are different 
for these two resonant states. Here we have shown that 
this is not the case: For PS doublets of single particle 
resonant states in nuclei, not only the energies, but also 
the widths are exactly the same in the PSS limit. 

Similar to what happens in bound states, when the 
PSS limit is not realized, the PSS in resonant states is 
broken. In the following we will use a solvable model to 
illustrate the conservation and the breaking of the PSS 
in resonant states. Note that this kind of investigation 
can also be done numerically with other potentials. 

We consider that S(r) and A(r) are both spherical 
square-well potentials, 



E(r) 
A(r) 



c, 


r <R, 


0, 


r>R, 


D. 


r<R, 


0. 


r > R, 



(9) 
(10) 



where C and D are constants and R is the width. For 
such potentials the wave function is continuous at r = R, 



but not its derivative because the derivative of a square- 
well potential is a delta-function. For r ^ R, Eq. (4) 
reads. 



dr 2 



1(1 + 1) 



- M+(r)M-(r) 



G(r) = 0. (11) 



The regular solution of this equation is just a combination 
of the Ricatti-Bessel functions, 



G(r) = (p/k) l+1 j l (kr), r<R 
i 
2 



G(r) = £ J°(p)h~{pr) J°{ P yhf(jpr) 



(12) 
,r> fl(13) 



with k = y/(e — C — M) (e-D + M). The coefficient 
(p/k) l+1 is inserted in accordance with Eq. (5) and 
ji(z) cx z l+1 as z — > 0. 

Next let us determine the Jost function from continu- 
ous conditions of radial wave functions at r = R, 



G(R+) = G{R-), F(R+)=F(R-), 



(14) 



where i?_ and R + mean that one approaches r = R 
from r < R and r > R respectively. A linear equation 
for the Jost function can be written immediately using 
the connection condition for G(r), 



i+i 



Ji(kR) = ~ J°(p)h~{pR) - J°{p)*h+(pR) 



(15) 

The continuous condition for F(r) can also be used to 
deduce a similar equation. The derivative of G(r) is not 
continuous at r = R, 



dG 
dr 



R4 



dG 
dr 



= -CF(R). 



(16) 



R- 



Approaching R from r < R, we can represent F(R) by 
1 



F(R) = 



M + C 

(p/k) I+1 
M + C-e 



\ dr r 



r=R 



kj'^kR) - -jj(kR)) . (17) 



By substituting Eqs. (12), (13), and (17) into Eq. (16) 
we get a linear equation for J^{p) and J^{p)* ■ Finally 
the Jost function for the small component reads, 



J?{P) 



P 



2ik l + : 
C 



ri (kR)phf\pR) - kj' r (kR)h±( P R) 
s — (kJi(kR) - ^Jl(kR)) hf(pR) 



(18) 

Now comparing Jost functions J^{p) and J7^(p), it is 
clear that they differ only in the part containing C be- 
cause they have the same I. In other words, in the PSS 
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FIG. 2. (Color online) The zeros of the Jost function 
on the complex energy plane in square-well potentials (9) and 
(10) with C = (solid symbols) and C = -66 MeV (half-filled 
symbols) for PS partners (diamond) andj'13/2 (square). 

The results with Woods-Saxon-like scalar and vector poten- 
tials are also shown as open symbols. 



limit, i.e., C = 0, we have J^{p) = <J$(p)- Conse- 
quently the PSS is conserved both in bound states and 
in resonant states. If C ^ 0, the PSS is broken and we 
can study the PS splitting of the energy and the width for 
resonant states. Due to the special form of the spherical 
square-well potentials, the PSS-breaking term is sepa- 
rated from the PSS-conserving term in the Jost function, 
which makes the study of the conservation or the break- 
ing of the PSS very convenient. 

The solution of J^{p) = can not be written in a 
closed form. Here because the Jost function is analytic 
near its zeros, the secant method can be used for search- 
ing the roots. Starting from an initial guess for a root, 
the iteration converges after a few steps. In Fig. 2 we 
show solutions in the complex energy plane for PSS dou- 
blets with I = 6, i.e., ft.xi/2 with k = —6 and ,713/2 with 
k' = 7 for square-well potentials with D = 650 MeV and 
R = 7 fm. In the PSS limit, i.e., C = 0, all the roots 
locate in the lower half plane and there are no bound 
states. We show in Fig. 2 three pairs of PS resonant 
partners by full diamonds and squares. The conservation 
of the PSS for single particle resonant states is clearly 
seen. When C = —66 MeV, there is one bound state 
only for hn/ 2 - Three pairs of PS partners of resonant 
states are shown by half-filled diamonds and squares. 
One finds the breaking of the PSS both in the bound 
states and in the resonant states. For PS doublets with 
other values of /, we observed similar behaviors concern- 
ing the exact conservation and the breaking of the PSS. 
We also studied resonances in Woods-Saxon-like poten- 
tials, W(r) = W /(l + exp[(r - R)/a\) (W = V or S) 
with parameters connected with 208 Pb given in Ref. [35]: 
the depths V -S = 650 MeV and V + S = -66 MeV, 
the diffusivity parameter a = 0.6 fm, and R = 7 fm. 
Resonance parameters are obtained with the real stabi- 
lization method [28]. The results are shown as open dia- 
monds and squares for hn/2 and J13/2, respectively. It is 



found that splittings of energy and width both become 
smaller compared with the results with the square-well 
potentials. The reason is that the derivative of S(r) is 
smaller due to a non-zero diffusivity parameter. 

Now we briefly discuss about protons. Due to the re- 
pulsive Coulomb interaction, S(r) can not be zero, nor 
its first derivative. Therefore the PSS limit can never be 
realized for protons. However, the Coulomb potential is 
relatively small compared with E(r) and the breaking of 
the PSS from the Coulomb interaction should be small. 

In summary, we show that the PSS in single particle 
bound and resonant states in nuclei can be investigated 
on the same footing within the relativistic framework by 
examining the zeros of Jost functions corresponding to 
small components of nucleon Dirac wave functions. In 
the PSS limit, i.e., the attractive scalar and repulsive 
vector potentials have the same magnitude but opposite 
sign, small components of PS doublets are exactly the 
same. Thus Jost functions describing the asymptotic be- 
havior of the radial wave functions arc identical to each 
other. When analytically continued to complex momen- 
tum plane, the resonant states, showing themselves as 
zeros of the Jost functions, are always paired in the PSS 
limit, which leads to the exact PSS. The conservation of 
the PSS in the PSS limit is also justified by examining 
the phase shift of continuum states. When leaving the 
PSS limit, the PSS in resonant states is broken. These 
conclusions arc tested for single particle resonances in 
spherical square-well and Woods-Saxon potentials. 
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